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We  describe  here  a  detailed  model  of  the  geometrical,  mate¬ 
rial,  and  electrical  properties  of  the  SCATHA  satellite  for  use 
with  the  NASA  Charging  Analyzer  Program  (NASCAP) .  Charging  cal¬ 
culations  in  an  intense  magnetospheric  substorm  environment  demon¬ 
strate  that:  (1)  long  booms  can  significantly  perturb  the  poten¬ 
tials  near  the  spacecraft,  and  (2)  discharging  by  sunlight  or  by 
active  control  can  cause  serious  time-dependent  differential 
charging  problems. 


1 .  INTRODUCTION 


We  have  developed  a  detailed  model  of  the  SCATHA  satellite 
for  use  with  the  NASA  Charging  Analyzer  Program  (NASCAP)  (refs. 

1  and  2) .  The  model  accounts  for  such  geometrical  complexities 
as  booms,  shadowing,  and  the  presence  of  insulating  materials 
over  portions  of  the  conducting  ground  of  the  space  vehicle.  The 
effects  of  photoemission  and  secondary  emission  caused  by  electron 
and  ion  impact,  active  control  devices  such  as  electron  and  ion 
beams,  and  surface  and  bulk  conductivity  are  included  in  the 
model.  To  our  knowledge,  this  model  represents  the  most  complete 
and  realistic  treatment  of  spacecraft  charging  attempted  to  date 
for  any  satellite. 

Section  2  below  describes  the  SCATHA  model  employed  in 
NASCAP.  A  detailed  shadowing  study  was  performed  for  a  geometri¬ 
cally  more  accurate  SCATHA  model;  this  work  is  described  in  Sec¬ 
tion  3.  We  have  performed  charging  calculations  for  one  environ¬ 
ment  using  the  present  model,  and  the  results  of  these  calcula¬ 
tions  are  described  in  Section  4.  Preliminary  conclusions  of  this 
study  are  summarized  in  Section  5. 


This  work  supported  by  the  National  Aeronautics  and  Space  Ad¬ 
ministration,  Lewis  Research  Center,  under  Contract  NAS3-21050. 
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2.  SCATHA  MODEL  DEVELOPMENT 


The  NASCAP  program  allows  the  specification  of  the  geometri¬ 
cal,  material,  and  electrical  properties  of  a  spacecraft  in  con¬ 
siderable  detail.  We  have  attempted  to  incorporate  the  most  cur¬ 
rent  and  complete  information  available  for  SCATHA  into  our  model. 
However,  the  present  model  is  meant  primarily  to  illustrate  the 
intended  level  and  scope  of  our  study,  rather  than  to  provide  the 
final  word  on  a  model  specification.  The  NASCAP  code  allows  model 
features  to  be  easily  altered  to  make  our  model  a  more  faithful 
representation  of  the  SCATHA  satellite  if  the  need  arises. 

Perspective  views  of  our  gridded  model  are  shown  in  figures  1 
and  2.  The  main  body  of  the  satellite  is  represented  as  a  right 
octagonal  cylinder,  with  the  aft  cavity  visible  in  figure  2.  The 
OMNI  antenna  and  the  SC9  cluster  of  experiments  are  visible  on 
the  forward  surface  of  the  satellite.  Our  model  reproduces  the 
actual  SCATHA  geometrical  features  extremely  well,  as  shown  in 
table  1.  Note  in  particular  that  the  treatment  of  booms  in 
NASCAP  allows  the  actual  boom  radii  to  be  reproduced  exactly  in 
the  model.  The  requirements  in  NASCAP  that  booms  parallel  coordi¬ 
nate  axes  and  intercept  mesh  points  in  all  grids  effectively  force 
any  long  booms  to  pass  through  the  center  of  the  innermost  mesh. 
Therefore,  our  present  model  includes  only  the  SC6 ,  SC11,  and  the 
two  SC2  booms,  with  the  orientations  fixed  at  right  angles  to  one 
another. 

Figure  3  illustrates  the  computational  space  in  which  NASCAP 
solves  Poisson's  equation  for  this  model.  Monopole  boundary  con¬ 
ditions  are  imposed  on  the  edges  of  the  outermost  grid,  which  is 
a  rectangular  prism  of  dimensions  1.6  *  1.6  x  3.2  m.  The  zone 
size  decreases  by  a  factor  of  2  in  each  of  the  four  successive 
inner  grids,  so  that  the  effective  resolution  is  11.5  cm  near  the 
satellite  body.  (Local  mesh  refinement  techniques  in  NASCAP  allow 
a  resolution  of  2.5  cm  for  selected  zones  on  the  satellite.) 

Our  model  includes  the  specification  of  15  distinct  exposed 
surface  materials,  each  of  which  is  specified  by  the  values  of 
some  13  user-supplied  parameters.  The  surface  materials  are  de¬ 
scribed  in  table  2.  We  have  attempted  to  find  experimentally 
measured  values  for  all  parameters;  where  this  has  not  been  pos¬ 
sible,  suitable  estimates  based  on  the  properties  of  similar  ma¬ 
terials  have  been  used.  Table  3  lists  the  values  employed  in  the 
calculations  reported  here.  The  analytical  expressions  in  which 
these  parameters  are  used  to  evaluate  net  surface  currents  are 
described  in  detail  in  reference  5.  The  formulation  of  electron 
backscattering  in  NASCAP  has  been  somewhat  modified  recently,  and 
the  newer  treatment  is  described  in  appendix  A.  The  exposed 
materials  are  illustrated  in  figure  4  in  which  the  locations  of 
several  of  the  SCATHA  experiments  are  also  shown.  Experiments  at 
the  ends  of  SCATHA  booms  are  modeled  as  a  single  boom  segment 
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whose  radius  is  adjusted  to  match  the  exposed  surface  area  of  the 
actual  experiment. 

The  model  includes  six  distinct  underlying  conductors: 
spacecraft  ground,  the  reference  band,  and  the  four  experiments 
SC2-1,  SC2-2,  SC6-1  and  SC6-2.  Each  of  these  underlying  con¬ 
ductors  is  capacitively  coupled  to  spacecraft  ground,  and  each 
can  be  separately  biased  with  respect  to  ground.  A  seventh  con¬ 
ductor  could  be  introduced  to  underlay  the  solar  cells  at  an  ap¬ 
propriate  bias.  In  this  study  the  reference  band  was  allowed  to 
float  and  all  other  conductors  were  biased  to  the  ground  potential. 

NASCAP  has  extensive  capabilities  to  model  particle  emitters 
and  detectors  located  on  the  spacecraft  body,  as  described  pre¬ 
viously  (ref.  2) .  These  features  of  NASCAP  can  be  used  in  the 
analysis  of  the  operation  of,  for  example,  the  SCATHA  experiments 
SC4,  SC5,  SC6,  SC7,  and  SC9 .  Such  studies  should  be  particularly 
helpful  in  determining  the  influence  of  spacecraft  fields  on 
particles  emitted  during  active  control,  and  in  determining  the 
source  of  particles  seen  at  detector  sites. 


3 .  SHADOWING  STUDY 


For  the  SCATHA  shadowing  study,  we  were  required  to  generate 
percent  shadowing  tables  for  various  experiments.  We  were  able 
to  generate  accurate  tables  using  relatively  small  amounts  of  com¬ 
puter  time:  less  than  5  minutes  Univac  1100/81  time  was  required 
for  a  table  of  7560  entries. 

Since  the  geometrical  capabilities  of  the  NASCAP  shadowing 
routines  are  more  general  than  the  rest  of  the  code,  we  were  able 
to  employ  a  SCATHA  model  for  shadowing  in  which  each  experiment 
was  treated  geometrically  in  much  finer  detail  than  in  the  model 
described  in  Section  2.  Figure  5  shows  the  level  of  detail  in  a 
perspective  view  of  the  ML12-7  experiments  on  the  forward  surface. 
Booms  were  placed  at  their  actual  locations  on  the  satellite,  and 
the  experiments  at  the  boom  ends  were  given  a  great  deal  of  geo¬ 
metrical  complexity.  Figure  6  shows  the  SC2-1,  SC1-4,  and  SC6-1 
booms  as  they  were  resolved  in  the  shadowing  study. 

These  detailed  geometrical  shapes  were  input  to  the  usual 
NASCAP  shadowing  routines  (HIDCEL)  for  table  generation.  The 
tables  cover  satellite  rotation  in  1°  increments  for  the  satel¬ 
lite  plane  deviations  from  the  sun  line  of  -5°  to  +5°. 


4 .  CHARGING  CALCULATIONS 


The  model  was  subjected  to  an  extremely  intense  substorm  de¬ 
scribed  by  a  superposition  of  two  Maxwellian  plasmas  with  the  fol¬ 
lowing  parameters: 
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The  effects  of  ambient  space  charge  were  neglected  in  the  solution 
of  Poisson's  equation  here,  since  the  mean  satellite  radius,  rg, 
is  much  smaller  than  the  plasma  Debye  length,  XD> 


r  ^  100  cm 
s  _ 

X„  7.  700  vl—  7,  2200  cm 


rs/XD  7,  0.05 


There  was  no  sunlight  present  in  the  first  calculation  described 
below. 


Potential  contours  during  the  initial  overall  charging  phase 
(^10“3  seconds)  are  shown  in  figures  7  and  8.  The  question  of 
whether  booms  have  a  significant  effect  on  the  sheath  potentials 
is  clearly  answered  by  examining  figure  9,  which  shows  potential 
contours  in  a  plane  a  half  meter  below  the  plane  of  the  booms. 
Figure  10  shows  similar  contours  in  a  calculation  with  the  booms 
omitted;  the  distortion  of  contours  by  the  booms  is  obvious. 

While  the  boom  radii  are  small,  v2  cm,  the  effect  on  potentials 
is  related  to  the  boom  capacitance,  which  varies  only  logarithmi¬ 
cally  with  radius.  This  results  in  long  range  potential  inter¬ 
actions  from  thin  booms,  where  the  characteristic  decay  distance 
is  closer  to  the  boom  length  than  to  the  boom  radius. 


The  rapid  initial  charging  is  followed  by  a  much  slower  de¬ 
velopment  of  differential  charging,  as  illustrated  in  figure  11. 
For  this  example  the  maximum  differential  developed  after  22  sec¬ 
onds  was  700  volts  and  the  maximum  field  strength  in  a  dielectric 
layer  was  24,000  volts/cm.  Figure  12  shows  contours  in  the  plane 
of  the  booms  after  22  seconds;  note  the  differential  charging  de¬ 
veloped  at  the  boom  ends  due  to  variations  in  the  material  proper¬ 
ties  between  the  experiments  and  the  boom  coatings . 
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The  two  Maxwellian  description  of  the  plasma  leads  to  a  low 
overall  charging  voltage  of  only  -7.3  keV  despite  the  presence  of 
a  plasma  component  with  an  electron  temperature  of  40  keV.  For 
the  particular  case  we  have  studied  here,  low  energy  protons  are 
being  collected  at  an  enormous  rate  and  these,  augmented  by  the 
secondary  electrons  they  produce,  balance  the  incident  electron 
current.  NASCAP  uses  a  proton  collection  model  in  which  the 
collection  increases  linearly  with  voltage,  which  is  valid  in  the 
present  case  where  rs/XD  is  small,  as  discussed  by  Laframboise 
(ref.  4) .  Table  4  shows  the  detailed  current  balance  near  equi¬ 
librium  for  the  boom  surface  material  in  the  presence  of  the 
double  Maxwellian  environment  described  above.  Also  shown  in 
table  4  is  a  similar  breakdown  for  the  same  material  subjected 
only  to  the  high  energy  single  Maxwellian  component.  The  equi¬ 
librium  potential  is  -32  keV  in  this  case,  indicating  that  the 
final  potentials  reached  would  have  been  much  lower  had  we  em¬ 
ployed  a  single  Maxwellian  plasma  model.  For  both  plasma  models, 
the  final  potentials  reached  will  depend  on  the  exact  values  em¬ 
ployed  for  the  proton  and  electron  induced  secondary  yields. 

Great  care  should  be  exercised  in  the  determination  of  the  values 
and  associated  error  estimates  for  parameters  which  affect  the 
production  of  secondary  electrons  in  these  and  similar  calcula¬ 
tions  . 

Finally,  the  atomic  number  dependence  of  backscatter  coeffi¬ 
cients  tends  to  make  high-Z  materials  charge  less  negatively  than 
other  elements.  For  SCATHA,  this  means  that  the  magnitude  of 
the  boom  potentials  will  be  significantly  lower  than  most  other 
surfaces,  since  exposed  platinum  constitutes  much  of  its  surface 
area. 


We  have  performed  a  similar  calculation  on  this  model  in 
which  the  sunlight  was  turned  on  after  22  seconds  of  charging  in 
eclipse.  The  photoemission  results  in  strong  differential  charg¬ 
ing  (^3  keV)  along  the  booms,  as  shown  in  figure  13.  In  our  model 
the  boom  surfaces  are  very  weakly  capacitively  coupled  to  the 
grounded  cable  shields  which  extend  the  length  of  the  booms,  while 
the  experiments  at  the  ends  of  the  SC2  and  SC6  booms  are  coupled 
closely  to  spacecraft  ground.  This  weak  coupling  has  the  effect 
of  allowing  the  booms  to  react  rapidly  to  environmental  perturba¬ 
tions  compared  to  the  rest  of  the  satellite,  leading  to  temporary 
conditions  of  high  differential  charging.  We  have  observed  simi¬ 
lar  effects  when  discharging  the  satellite  with  an  electron  gun. 

The  potentials  near  the  satellite  in  sunlight  are  dominated 
by  the  monopole  field  of  the  spacecraft  body.  A  photoemitting 
boom  surface  element  can  discharge  only  to  the  value  of  the  local 
monopole  potential,  since  further  discharge  is  limited  by  immedi¬ 
ate  reflection  of  photoelectrons.  This  has  the  amazing  conse¬ 
quence  that  the  booms,  strongly  perturbing  in  eclipse,  now  seem 
to  disappear  in  the  potential  contours  near  the  satellite  body. 
Note  that  significant  differential  charging  in  sunlight  along  the 


SC2  booms  will  certainly  persist  at  equilibrium  due  to  large  dif¬ 
ferences  between  the  photoemission  from  surfaces  on  booms  and  on 
the  SC2-1  and  SC2-2  experiments.  Our  calculations  neglect  any 
effective  surface  conductivity  parallel  to  the  booms  due  to  the 
presence  of  a  photosheath.  The  surface  conductivity  features  of 
NASCAP  could  easily  be  invoked  to  simulate  this  effect,  which 
would  reduce  the  magnitude  of  the  differential  charging  observed 
here. 


The  calculations  reported  here  were  performed  on  the  Univac 
1100/81  computer  at  Systems,  Science  and  Software.  Each  cycle  of 
charging  and  solution  of  the  potential  equations  required  approxi¬ 
mately  15  minutes  CPU  time  during  differential  charging,  and  5 
minutes  CPU  time  when  no  differential  charging  occurred.  Approxi¬ 
mately  10  cycles  of  each  type  were  required  for  the  calculations 
reported  here.  We  have  developed  a  second  SCATHA  model  for  test¬ 
ing  purposes  in  which  the  zone  size  is  twice  that  of  the  model 
presented  here  and  the  booms  are  shortened;  computer  times  are 
reduced  by  roughly  80  percent  for  this  model,  and  all  of  the  re¬ 
sults  described  above  can  be  observed  in  calculations  using  the 
smaller  model.  The  half-scale  model  will  be  useful  whenever  fine 
resolution  on  the  satellite  surfaces  is  not  required. 


5.  CONCLUSIONS 


We  have  completed  the  development  of  a  detailed  model  of  the 
SCATHA  satellite.  Preliminary  results  from  calculations  in  one 
magnetospheric  environment  indicate  that: 

•  The  presence  of  a  low  energy  component  in  a  two  Maxwellian 
description  of  the  magnetospheric  environment  reduces  the 
maximum  charging  of  a  satellite  relative  to  that  found  for 
a  single  Maxwellian. 

•  The  booms  have  substantial  impact  on  potentials  near  the 
spacecraft  in  eclipse. 

•  The  use  of  high  atomic  number  coatings,  such  as  platinum 
on  the  booms,  may  increase  the  severity  of  differential 
charging. 

•  Discharging  by  sunlight  or  by  active  control  may  lead  to 
transient:  increases  in  differential  charging  along  the 
booms  due  to  the  weak  coupling  of  the  booms  to  spacecraft 
ground. 

Our  calculations  demonstrate  that  the  prediction  of  space¬ 
craft  potentials  for  SCATHA  is  an  exceedingly  complex  problem,  in 
which  the  full  capabilities  of  the  NASCAP  treatment  of  geometrical 
features,  material  properties,  and  dynamic  interaction  with  the 
environment  are  utilized.  We  plan  to  continue  this  study  of 
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SCATHA  using  NASCAP  with  particular  emphasis  on  boom  perturbations 
and  the  effects  of  active  control. 


APPENDIX  A.  ELECTRON  BACKSCATTER 


Electron  backscatter  is  modeled  in  NASCAP  as  a  function  of 
electron  energy  and  mean  atomic  number  of  backscattering  material. 
The  formulation  first  used  in  NASCAP  (ref.  5)  was  valid  only  for 
low  Z-materials.  To  remove  this  restriction  we  have  used  a  for¬ 
mula  of  Burke  (ref.  6)  to  obtain  the  backscatter  coefficient  for 
isotropically  incident  electrons  as 

nx  =  0.475  Z0*177  -  0.40  .  (Al) 

The  backscatter  coefficient  for  normal  incidence,  n  ,  is  then 
found  by  solving  the  equation 

nl  =  2[1  ‘  r'o(1'ln  nQ)]/(ln  n0)2  ( A2 ) 

which  comes  from  assuming  the  angular  dependent  backscatter  co¬ 
efficient  (ref.  7)  to  be 

n(0)  =  nQ  exp [- ( £n  nQ)  (1  -  cos6)]  . 

The  energy  dependence  (ref.  4)  is  then  taken  to  be 

nQ(e)  =  YU)  (nQ  +  0.1  exp(-e/5)) 

!0  e  <  50  eV 

Jin  (20  e)/Jln  20  50  eV  <  e  <  1  keV 

1  e  >  1  keV 

where  e  is  in  keV. 

The  energy  dependent  nQ  from  (A4)  is  then  used  in  (A2)  or 
(A3)  to  calculate  the  relevant  backscatter  coefficient. 


(A3) 

(A4) 
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TABLE  1.  COMPARISON  OF  ACTUAL  SCATHA  GEOMETRICAL  FEATURES  TO 
GRIDDED  NASCAP  MODEL 


Zone  Size  =  4.54  in.  (11.5  cm) 


SCATHA 

MODEL 

Radius 

33.6  inches 

32.0  inches 

Height 

68.7 

68.0 

Solar  Array  Height 

29 

27.2 

Bellyband  Height 

11.3 

13.6 

SC9-1  Experiment 

9.2  x  6  x  8 

9.1  x  4.5  x  9.1 

SC6-1  Boom 

1.7  (radius) 

1.7 

118  (length) 

113.2 

Surface  Area 

4 

2.16  x  io  sq.  in. 

4 

2.11  x  10  sq. 

Solar  Array  Area 

1.23  x  io4 

1.15  x  104 

Forward  Surface  Area 

0.36  x  104 

0.34  x  104 
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TABLE  2.  EXPOSED  SURFACE  MATERIALS 


GOLD: 

SOLAR: 

WHITEN: 

SCREEN: 

YELOWC : 
GOLDPD : 

BLACKC : 
KAPTON : 
SI02: 
TEFLON : 
INDOX: 
YGOLDC : 

ML12: 

ALUM: 
BOOMAT : 


gold  plate 

solar  cells,  coated  fused  silica 

non-conducting  white  paint  (STM  K792) 

SC5  screen  material,  a  conducting  fictitious 
material  which  absorbs  but  does  not  emit 
charged  particles 

conducting  yellow  paint 

88  percent  gold  plate  with  12  percent  conductive 
black  paint  (STM  K748)  ’ n  a  polka  dot  pattern 

conductive  black  paint  (STM  K748) 

kapton 

Si02  fabric 

teflon 

indium  oxide 

conducting  yellow  paint  (50  percent) 
gold  (50  percent) 

ML12-3  and  ML12-4  surface,  a  fictitious  material 
whose  properties  are  an  average  of  the  proper¬ 
ties  of  the  several  materials  on  the  ML12  sur¬ 
faces 

aluminum  plate 
platinum  banded  kapton 
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TABLE  3.  MATERIAL  PROPERTIES  FOR  EXPOSED  SURFACES 
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TABLE  3.  (Continued) 

aThe  materials  are  described  in  Table  2. 

^The  thirteen  properties  are  as  follows  (see  Reference  4  and 
Appendix  A  for  further  details) : 

Property  1:  Relative  dielectric  constant  for  insula¬ 

tors  (dimensionless) . 

Property  2:  Thickness  of  dielectric  film  or  vacuum  gap 

(meters) . 

Property  3:  Electrical  conductivity  (mho/m).  The  value 

«  indicates  a  vacuum  gap  over  a  conducting 
surface. 

Property  4:  Atomic  number  (dimensionless). 

Property  5:  Maximum  secondary  electron  yield  for  elec¬ 

tron  impact  at  normal  incidence  (dimen¬ 
sionless)  . 

Property  6 :  Primary  electron  energy  to  produce  maximum 

yield  at  normal  incidence  (keV) . 

Properties  7-10:  Range  for  incident  electrons.  Either : 

P  P 

8  10 
Range  =  P?E  °  +  PgE 

where  the  range  is  in  angstroms  and  for 
the  energy  in  keV, 

or 

P-7  =  -1.  to  indicate  use  of  an  empirical 
range  formula 

Pg  =  density  (g/cm3) 

P10  =  mean  atomic  weight  (dimensionless) . 

Property  11:  Secondary  electron  yield  for  normally 

incident  1  keV  protons. 

Property  12 :  Proton  energy  to  produce  maximum  secondary 

electron  yield  (keV) . 

Property  13:  Photoelectron  yield  for  normally  incident 

sunlight  (A/m2) . 

The  dielectric  constant  and  thickness  for  the  boom  surfaces  were 
chosen  to  reflect  the  effective  capacitance  to  the  underlying 
cable  shield. 
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TABLE  4.  COMPONENTS  OF  INCIDENT  AND  EMITTED  CURRENTS  (10~5  A/m2) 
FOR  BOOM  SURFACE  MATERIAL  NEAR  STEADY  STATE. 


Double 

Maxwellian 

Single 

Maxwellian 

Potential 

-7000  Volts 

-32,000  Volts 

Incident  Electrons 

-4.6 

-2.3 

Resulting  Backscatter 

2.7 

1.4 

Resulting  Secondaries 

.7 

.4 

Incident  Protons 

.6 

.2 

Resulting  Secondaries 

.6 

.  3 

Figure  3.  Computational  space  surrounding  the  SCATHA  model,  show¬ 
ing  the  nesting  of  the  grids.  The  tic  marks  along  the  axes  indi¬ 
cate  the  outer  grid  zone  size;  the  zone  size  decreases  by  a  factor 
of  two  in  successive  grids. 
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Figure  4a.  SCATHA  model  with 
exposed  surface  materials  il¬ 
lustrated. 
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Figure  4b.  SCATHA  model  with 
exposed  surface  materials  il¬ 
lustrated  . 
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Figure  4c.  SCATHA  model  with 
exposed  surface  materials  il¬ 
lustrated. 


Figure  7 .  Potential  contours  in  a  vertical  plane  through  SCATHA 
center  (only  two  of  the  four  grids  are  plotted) .  Note  the  con¬ 
tours  extending  into  the  aft  cavity.  Time  ^10~3  seconds.  Con¬ 
tours  from  -450  to  -1250  volts  in  50  volt  steps. 
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Figure  8.  Potential  contours  in  a  horizontal  plane  through  SCATHA 
center.  Time  ^10"3  seconds.  Contours  from  -300  to  -1200  volts  in 
100  volt  steps.  The  relative  orientations  of  the  booms  is  the 
same  in  later  figures.  The  dimples  in  the  potential  contours 
near  the  boom  ends  are  artifacts  associated  with  an  imperfect 
match  of  potential  interpolation  functions. 
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'  X 

Figure  9.  Potential  contours  in  a  horizontal  plane  1  m  below 
SCATHA  center.  Time  ^10-3  seconds.  Contours  from  -250  to  -1150 
volts  in  50  volt  steps. 


Figure  10.  Potential  contours  in  a  horizontal  plane  1  m  below 
SCATHA  center  for  a  model  in  which  the  booms  have  been  removed. 
Time  ^10-3  seconds.  Contours  from  -300  to  -1900  volts  in  100 
volt  steps. 


TIME  (SEC) 

Figure  11.  Spacecraft  potential  versus  time  for  two  points  on 
SCATHA  satellite. 


Figure  12.  Potential  contours  in  a  horizontal  plane  through 
SCATHA  center,  with  differential  charging  along  booms.  Time 
^22  seconds.  Contours  from  -2000  to  -7000  volts  in  500  volt 
steps . 


Figure  13.  Potential  contours  for  sunlit  case  in  a  horizontal 
plane  through  SCATHA  center.  Time  ^38  seconds.  Contours  from 
-1000  to  -7500  volts  in  500  volt  steps. 


